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Abstract 

We study the asymptotic behavior of wavelet coefficients of random processes with long 
memory. These processes may be stationary or not and are obtained as the output of 
non-linear filter with Gaussian input. The wavelet coefficients that appear in the limit are 
random, typically non-Gaussian and belong to a Wiener chaos. They can be interpreted 
as wavelet coefficients of a generalized self-similar process. 
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1. Introduction 

Let X = {Xn}n& be a stationary Gaussian process with mean zero, unit variance and 
spectral density /(A), A G (— vr, vr] and thus covariance equal to 



r{n) = E(XoX„) = T e^"V(A)dA 

J — 7r 



The process {Xn}n& is said to have short memory or short-range dependence if /(A) is 
bounded around A = and long memory oi long-range dependence if /(A) — t- oo as A — )■ 0. 
We will suppose that {Xn}n& has long-memory with memory parameter d > 0, that is, 

/(A) ~ \X\-^^f*{X) as A ^ 
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where /*(A) is a bounded spectral density which is continuous and positive at the origin. 
It is convenient to interpret this behavior as the result of a fractional integrating operation, 
whose transfer function reads A i— )■ (1 — e~''^)~'^. Hence we set 

f{\) = \l-e--^'\-'''n\), \ei-n,n]. (1) 

We relax the above assumptions in two ways : 

1. Consider, instead of the Gaussian process {Xn\n& the non-Gaussian process 

where G is a non-linear filter such that E[G(X„)] = and E[G(X„)^] < oo. The 
non-hnear process {G{Xn)}nei is said to be subordinated to the Gaussian process 

2. Drop the stationarity assumption by considering a process {Yn}nez which becomes 
stationary when differenced i^' > times. 

We shall thus consider {Yn}n& such that 

where (AF)„ = F„ — and where {Xn}n& is Gaussian with spectral density / satisfy- 
ing (HD. 

Since Y = {Yn}nei is random so will be its wavelet coefficients {Wj^k, j ^ 0, k & 
which are defined below. Our goal is to find the distribution of the wavelet coefficients at 
large scales j — t- oo. This is an important step in developing methods for estimating the 
underlying long memory parameter d. The large scale behavior of the wavelet coefficients 
was studied in [H in the case where there was no filter G, that is, when F is a Gaussian 
process such that A^Y = X, and also in the case where F is a non-Gaussian linear process 
(see 3). 

We obtain our random wavelet coefficients by using more general linear filters that 
those related to multiresolution analysis (MRA) (see for e.g. [sj, Q). In practice, however, 
the methods are best implemented using Mallat's algorithm and a MRA. Our filters are 
denoted hj where j is the scale and we use a scaling factor ■jj ^ oo as j ^ oo. In the case of 
a MRA, 7j = 2^ and hj are generated by a (low pass) scaling filter and its corresponding 
quadratic (high pass) mirror filter. More generally one can use a scaling function (p and a 
mother wavelet ip to generate the random wavelet coefficients by setting 



dt, (2) 



where V'j.fc = 2-^/2^(2-^t - A;), j > 0. Observe that we use here the engineering convention 
that large values of j correspond to large scales and hence low frequencies. If (f and ip have 
compact support then the corresponding filters hj have finite support of size 0{2^). For 
more details on related conditions on Lp and ip (see [H). 
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The idea of using wavelets to estimate the long memory coefficient d goes back to 
Wornell and al. (fl) and Flandrin (0, 0, B 8). See also Abry and al. ([lo|,lll|). Those 
methods are an alternative to the Fourier methods developed by Fox and Taqqu ((l2j) 



and Robinson 
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14| . For a general comparison of Fourier and wavelet approach, see 
The case of the Rosenblatt process, which is the Hermite process of order q = 2, was 



studied by [16 



The paper is structured as follows. In Section |2l we introduce the wavelet filters. The 
processes are defined in Section [3] using integral representations and Section H] presents the 
so-called Wiener chaos decomposition. The main result and its interpretations is given 
in Section [5l It is proved in Section El Auxiliary lemmas are presented and proved in 
Sections [7] and M 



2. Assumptions on the wavelet filter 

The wavelet transform of Y involves the application of a linear filter hj{r),T G Z, at 
each scale j > 0. We shall characterize the filters hj by their discrete Fourier transform : 

/^.(A) = $^/i,(r)e~^^^AG[-7r,7r]. 
Assumptions on hj are stated below. The resulting wavelet coefficients are defined as 

where 7j ^ oo is a sequence of non-negative scale factors applied at scale j, for example 
7j = 2K We will assume that for any m G Z, 

lim^ = 7^>0. (3) 

As noted, in this paper, we do not assume that the wavelet coefficients are orthogonal 
nor that they are generated by a multiresolution analysis. Our assumptions on the filters 
hj are as follows : 

a. Finite support: For each j, {hj{T)}rez has finite support. 

b. Uniform smoothness : There exists M > a> 1/2 and C > such that for all j > 

and A G [— tt, tt], 

- (l + 7,|A|)A^+° ■ 
By 27r-periodicity of hj this inequality can be extended to A G M as 



(4) 



'^^^^^'-^(l + 7,|{A}|r^- 
where {A} denotes the element of (— vr, vr] such that A — {A} G 27rZ. 
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c. Asymptotic behavior: There exists some non identically zero function hoo such that 
for any A G M, 

.lim {i;'%{i;'x)) = h^{x) . (6) 



Observe that while hj is 27r-periodic, the function hoo is a non-periodic function on R 
(this follows from (|T2|) below). For the connection between these assumptions on hj and 
corresponding assumptions on the scaling function ip and the mother wavelet ip in the 

classical wavelet setting ([2]) (see [1]). In particular, in that case, one has h^o = 'p{0)ip. 
Our goal is to study the large scale behavior of the random wavelet coefficients 

Wj,k = J2 h^i^^f" - ^)^^ = E ^^-(^^-^ - ^) {^'''GiX))^ , (7) 

where we set symbohcally Ye = {A-^G{X))^ for {A^Y)i = G{Xi). 

By Assumption (jl]), hj has null moments up to order M — 1, that is, for any m G 
{O,--- ,M-1}, 

J2h,ii)r = o. (8) 

Therefore, since M > K, hj can be expressed as 

%{X) = {l~e~'Yhf\X), (9) 
where h^j^^ is also a trigonometric polynomial of the form 

hf'W = j:hf\r)e-'\ (10) 

tGZ 

since hj has finite support for any j. Then we obtain another way of expressing Wj^k, 
namely, 

W,,k = J2hf\^,k-i)G{Xe). (11) 

^ez 

We have thus incorporated the linear filter in (jTj) into the filter hj and denoted the 
new filter h^j^\ 

Remarks 

1. Since {G(X^),£ G Z} is stationary, it follows from (fTTjl that {Wj^k, A; G Z} is station- 
ary for each scale j. 

2. Observe that A^Y is centered by definition. However, by ([8]), the definition of Wj^k 
only depends on A'^Y. In particular, provided that M > K + 1, its value is not 
modified if a constant is added to A^Y, whenever M > K + 1. 
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3. Assumptions (jl]) and imply that for any A G 



lA 



M 



IMA)|<C^^^^^. (12) 



Hence h^o e L^(M) since a > 1/2. 
4. The Fourier transform of /, 

mm= [ f{t)e-''"^dH, eeM^ (is) 

is defined for any / G L^(R'', C). We let h^o be the L^(M) function such that h^o = 

3. Integral representations 

It is convenient to use an integral representation in the spectral domain to represent the 



ipre 

random processes (see for example |l7l. Il8[|). The stationary Gaussian process {Xk, k E Z,} 



with spectral density ([T]) can be written as 

X,= e'^'f'\X)m{\)= I ^_ _il|; dW^(A), I en. (14) 

J —IT J —IT \^ ^ I 

This is a special case of 



I{g) = / g{x)m{x), (15) 

where W{-) is a complex-valued Gaussian random measure satisfying 

E(I?(A)) = for every Borel set A in M , (16) 

^{W{A)W{B)) = \AnB\ for every Borel sets A and 5 in M , (17) 

n n 

^ W{Aj) = l?(|J Aj) if Ai, ■ ■ ■ , An are disjoint Borel sets in R , (18) 



W{A) = W{-A) for every Borel set A in M . (19) 
The integral (fT5|) is defined for any function g G L^(R) and one has the isometry 

Ei\Tig)\')= [ \g{x)\Mx. 
Jr 

The integral I{g), moreover, is real-valued if 



9{x) = g{-x) 



We shall also consider multiple Ito- Wiener integrals 

where the double prime indicates that one does not integrate on hyperdiagonals Aj = 
±Aj, i 7^ j. The integrals Iq{g) are handy because we will be able to expand our non-linear 
functions G{Xk) introduced in Section [T] in multiple integrals of this type. 

These multiples integrals are defined as follows. Denote by L'^{M.'^,C) the space of 
complex valued functions defined on M*^ satisfying 

g{~xi, ■ ■ ■ , -Xg) = g{xi, ■ ■ ■ , Xg) for (xi, ■ ■ ■ , Xg) E M" , (20) 

llfi'lli^ / \g{xi,- ■ ■ ,Xg)f dxi- ■ -dxg < oo . (21) 

Let L^{W,C) denote the set of functions in 1?{WX) that are symmetric in the sense 
that g = g where g{xi,--- ,Xq) = l/o'! ^^.^'(^^^{i), ■ " " ^^aiq)), where the sum is over all 
permutations of {l,...,g}. One defines now the multiple integral with respect to the 
spectral measure W hj a. density argument. For a step function of the form 

q = > Cj, ... 1 1 A X ■ ■ ■ X 1a 

je=±l,-,±N 

where the c's are real-valued, Aj^ = — A_j^ and Aj^ fl Aj^ = if £ 7^ m, one sets 

%i9) = T^'^'nr-.MWiA,,) . ■■WiA.J . (22) 

je=±l,-,±N 

Here, indicates that one does not sum over the hyperdiagonals, that is, when ji = ±jm 
for £ m. The integral Ig verifies that 

EaA90/,.te))={?'<fr;;l^;; (23) 

Observe, moreover, that for every step function g with q variables as above 

Since the set of step functions is dense in L'^{W, C), one can extend Iq to an isometry from 
L'^{M.'^, C) to L^(f2) and the above properties hold true for this extension. 

Remark. Property fl20l) of the function / in L'^{M.'^, C) together with Property f|T9l) of W 
ensure that Iq{f) is a real-valued random variable. 
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4. Wiener Chaos 

Our results are based on the expansion of the function G, introduced in Section [H in 
Hermite polynomials. The Hermite polynomials are 

2 rf" / .2 



2 (7^ / 
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in particular, Hq{x) = l,Hi{x) = x,H2{x) = x"^ — 1. If X is a normal random variable 
with mean and variance 1, then 



E{H,{X)H,iX)) = [ Hg{x)H,,{x)^e-'^'^^dx = q\S, 

Moreover, 



q,q ■ 



+ 00 

Cr, 



G{X) = J2^H,{X), (24) 

where the convergence is in L'^{Q) and where 

c, = E{G{X)H,{X)) . (25) 

The expansion (l24l) is called a Wiener chaos expansion with each term in the chaos expan- 
sion living in a different chaos. The expansion (!24l) starts at g = 1, since 

Co = E{G{X)Ho{X)) = E{G{X)) = , 

by assumption. The condition E{G{Xy) < oo implies 

+00 2 

E^<^- (26) 

Hermite polynomials are related to multiple integrals as follows : if X = J^g{x)dW{x) 
with E(X^) = J^\g{x)\'^dx = 1 and g{x) = g{—x) so that X has unit variance and is 
real-valued, then 

H,{X) = T,{g^'^) = [ g{xi) ■ ■ ■ g{x,)dW{xi) ■ ■ ■ dW{x,) . (27) 

The expansion of G induces a corresponding expansion of the wavelet coefficients 
Wj^k, namely, 

+ 00 

W.. 

^1 j,k ' 

q=l ^ 

where by (fTTj) one has 



+ 00 

q=l 



w\$ = Y,hf\l,k-E)H,{X,). (29) 



The Gaussian sequence {Xn}n& is long-range dependent because its spectrum at low 
frequencies behaves like lAI^^"^ with d > and hence explodes at A = 0. What about the 
processes {Hq{X£)}£ for q >27 What is the behavior of the spectrum at low frequencies? 
Does it explodes at A = 0? The answer depends on the respective values of q and d. Let 
us define 

gc = max{gGN : q<l/{l-2d)} , (30) 

and 

d{q) = qd+{l-q)/2. (31) 

One has 

d{q) > if g < gc, that is if g < 1/(1 - 2d) . (32) 

The following result shows that the spectral density of {Hq{Xi)}i^z has a different behavior 
at zero frequency depending on whether q < qc or q > qc- It is long-range dependent when 
q < qc and short-range dependent when q > q^. We first give a definition. 

Definition 4.1. The convolution of two locally integrable {2tt) -periodic functions gi and 
g2 is defined as 

/•TT 

(fi-i *fi'2)(A) = / gi{u)g2{X- u)du . (33) 



Moreover the q times self- convolution of g is denoted by g^*'^\ 

Lemma 4.1. Let q be a positive integer. The spectral density of {Hq{X£)}£^i is 

g!/('^«)=g!(/^...^/), 

where the spectral density f of {Xi}^^^ is given in (Qp. Moreover the following holds : 

(i) If Q ^ Qc, then X^^(i^ f(*i^(^X) is bounded on \ E (0,vr) and converges to a positive 
number as X 10. 

(a) If q > qc, then f^*'^\X) is bounded on X E (0,7r) and converges to a positive number 
as X 10. 

Hence if q < qc, {IIq{Xe)}e has long memory with parameter d{q) > whereas if q > qc, 
{Hq{X)}i has a short-memory behavior. 

Proof. By definition of Hq and since X has unit variance by assumption, we have 
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EiHq{Xe)HqiX,+J) = q\ /(A)e^^"^dA 



Using the fact that, for any two locally integrable (27r)-periodic functions gi and g2, one 
has 

'^^7i*^?2)(A)e^^"^dA= r g,{u)e'-'^du X f g^iv^'^'^dv , 

r J —TV J —TV 
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we obtain that the spectral density of {Hg(Xi)}e is q\f^*'^\ 

The properties of f^*'^^ stated in Lemma |4?T] are proved by induction on q using Lemma [8l2l 
Observe indeed that if Pi = d{q) and P2 = 2(i, then 

/3i + /32 - 1 = 2d{q) + 2d-l = {2dq + 1 - g) + 2rf - 1 = 2(g + - (g + 1) + 1 = 2d{q + 1) . 

□ 

Now, consider the expansion of A^Yg = G{Xe) = '}2^=qS^qll^-)Hq{Xf), where 

go = min{g > 1, c, ^ 0} . (34) 

The exponent go is called the Hermite rank of A^Y. 

In the following, we always assume that at least one summand of A^Yf has long memory, 
that is, in view of Lemma I4.H 

go < gc • (35) 



5. The result and its interpretations 

In this section we describe the limit in distribution of the wavelet coefficients { Wj+m,fc}m,fc 
as j — )■ cxD, adequately normalized, and we interpret the limit. Recall that Wj^rn,k involves 
a sum of chaoses of all order. In the limit, however, only the order go will prevail. The 
convergence of finite-dimensional distributions is denoted by 

Theorem 5.1. As j — )■ 00, we have 

^-(%o)+x)^^.^^^^^ m. A; G Z} ^1 (/^(O))^^/^ {yj.'^^/^ m. A: G z} , (36) 
where for every positive integer q, 

-^im) y^,(i(c, + ... + g)^ lCi|^---|C,l^ ^^'^ dW^lC). (37) 
This Theorem is proved in Section El 



Interpretation of the limit. 



The limit distribution can be interpreted as the wavelet coefficients of a generalized 
Hermite process defined below, based on the wavelet family 

{/ioo,m,fc(t) = 7;;'/'/ioo(-7™ + k), m,keZ} . (38) 

This wavelet family is the natural one to consider because the Fourier transform h^oiX) is 
the rescaled limit of the original hj{\) as indicated in (El). 

A generalized process is indexed not by time but by functions. The generalized Hermite 
processes for any order g in {1, . . . , gc} are defined as follows : 
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Definition 5.1. Let < d < 1/2 and let q be a positive integer such that < q < 1/(1— 2d) 
and K > 0. Define the set of functions 



0(0 ' |e|«-^-2''^-2^de < OO 



where 6 = The generalized random process Z^^^ is indexed by functions 9 G S^^'^ and 
is defined as 



where 6 = '^[9] as defined in l[T3\) . 

Now fix (m, k) G and choose a function hoo^rn,k{t),t G M as in fl38l) . so tliat 



d[hoo,m,km = d[i;n'^'hU--f;^H + kM) = iiJ'^' e--^^-^ h^i-f^O ■ (40) 

Lemma 5.1. The conditions on d and q in Definition \5.1\ ensures the existence of zl^^{9). 
In particular, 

hoc,m,k e 5^5^ for allK e{0,...,M} , 

and hence Z^^^ {hoo,m,k) is well-defined. 

Tfiis Lemma is proved in Section [71 

By setting in fl39|) . 9 = /?.oo,m,fc, defined in fHOl) . we obtain for all (m, k) G Z^, 

^m,k ~ ^q,d \"'oo,m,k) ■ 

Hence tlie riglit-liand side of (136|1 are tlie wavelet coefficients of tlie generalized process 
^g^'' with respect to the wavelet family {hoo,m,k, rn,k G Z}. 

In the special case q = 1 (Gaussian case), this result corresponds to that of Theo- 
rem 1(b) and Remark 5 in [ij, obtained in the case where jj = 2K In this special case, we 
have z'^^ = Bi^d+K) , where is the centered generalized Gaussian process such that for 
all 9,,92eS[% 

Cov(i?(,)(ei),i?(,)(^^2))= / \\\-'%{\)92{X} dX . 

It is interesting to observe that, under additional assumptions on 9, for K > 1, Z^^^{9) 
can also be defined by 

^i?w= / ^iy(t)^dt, (41) 

where {Z^^^{t), t G M} denotes a measurable continuous time process defined by 

Z^Vit) = / dW{u,) ■ ■■dW{u,), teR. (42) 



(i(ui H h Ug))^\ui ■ --u. 



■<?i 
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If, in (HTj) we set = 1, we recover the usual Hermite process as defined in 19| wfiich has 



stationary increments. The process Z^^^ (t) can be regarded as the Hermite process Z^^j (t) 
integrated K — 1 times. In the special case where i^' = g = 1, we recover the Fractional 
Brownian Motion {BH{t)}t&R with Hurst index H = d + 1/2 G (1/2, 1). 

In the case = we cannot define a random process zj^^j (t) as in (1421) . The case K = 
would correspond to the derivative of the Hermite process Z^^j (t) but the Hermite process 
is not differentiable and thus the process Z^^l{t),t G M is not defined. When K = one 

can only consider the generalized process Z^°j(6'). Relation fH2|) can be viewed as resulting 
from fl39l) and fHTj) by interverting formally the integral signs. 
We now state sufficient conditions on 6 for fj4T]) to hold. 

Lemma 5.2. Let q be a positive integer such that < g < 1/(1 — 2d) and K > 1. Suppose 
that 9 G Sg^^ is complex valued with at least K vanishing moments, that is, 

[ e{t) dt = for all £ = 0, 1, . . . , K - 1 . (43) 

Suppose moreover that 

f \e{t)\\t\^+^'^-'^^^i dt<oo . (44) 
Jr 

Then Relation [J^H holds. 

This lemma is proved in Section [7l 

If, for example, the hj are derived from a compactly supported mult iresolut ion analysis 
then hoo will have compact support and so hoo,m,k will satisfy ( Hil) . In this case, the limits 
Y^'^^ in Theorem 15 . 1 1 can therefore be interpreted, for m, G Z as the wavelet coefficients 

of the process Z^^^ belonging to the g-th chaos. This interpretation is a useful one even 
when the technical assumption (14^ is not satisfied. 

Self-similarity. 

The processes Z^^^ and zj^^^ are self-similar. Self-similarity can be defined for processes 
indexed by t G M as well as for generalized processes indexed by functions 6 belonging to 
some suitable space S, for example the space S^^^ defined above. 

A process {Z{t), t G M} is said to be self-similar with parameter if > if for any a > 0, 

{a^ Zitja), t G M} = {Z{t), t G M} , 

where the equality holds in the sense of finite-dimensional distributions. A generalized 
process {Z{6), 6 E S} is said to be self-similar with parameter if > if for any a > and 

ees, 

Z{9'''") = Z{6) , 
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where 6°"'^ {u) = a ^0{u/a) (see [l7|, Page 5). Here S is assumed to contain both O""'^ and 

e. 

Observe that the process {Z^^^{t)^ t G M}, with > 1 is self-similar with parameter 
H = K + qd-q/2 = {K-l) + {d{q) + 1/2) . (45) 

As noted above zj^^^ can be regarded as Z^^j integrated K — 1 times. 

The generalized process {zjj^\9), 6 G Sj^^^}, which is defined in (1391) with > 0, is 
self-similar with the same value of H as in (BSl) . but this time the formula is also valid for 
K = 0. 

In particular, the Hermite process {K = 1) is self-similar with H = d{q) + 1/2 G (1/2, 1) 
and the generalized process Zg°j(^^) with i^' = is self-similar with H = d{q) — 1/2 G 
(-1/2,0). 

Interpretation of the result. 

In view of the preceding discussion, the wavelet coefficients of the subordinated process 
Y behave at large scales (7^ — )■ 00) as those of a self-similar process zj^^^ living in the 
chaos of order go (the Hermite rank of G) and with self-similar parameter K + (i(go) ~ 1/2. 

6. Proof of Theorem 15.11 

Notation. It will be convenient to use the following notation. We denote by Eg, q > 1, 
the C — )■ C function defined, for aA\ y = {yi, . . . , yg) by 



^M = Y.y^- (46) 



i=l 



With this notation Y^'^^ in Theorem 15.11 can be expressed as 

where o denotes the composition of functions. 

We will separate the Wiener chaos expansion ( !28|) of Wj^k into two terms depending on 
the position of q with respect to qc- The first term includes only the g's for which Hg{x) 
exhibits long-range dependence (LD), that is, 

q=0 ^' 
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and the second term includes the terms which exhibit short-range dependence (SD) 

oo 

<?=9c + l ^' 

Using Representation (fT4l) and ( |27ll since X has unit variance, one has for any £ G Z, 
= I expoS,(i£0 X (r^(0)'^' dl?(ei)---dl?(e,) . 

J (-7r,7r]'3 

Then by (I29]),([l0]) and (jS]), we have 



J2 hf\l,k -i) [ exp oE,(i£0 X (/«^(0)'^' dl?(ei) ■ ■ ■ dW{Q 
/ f$^/^r(7,A:-^)expoS,(i£0 1 (^^(0)'^' dl?(a) ■ ■ ■ dW?(e,) 
/■ e^'^'^^'^^) (5^/if)(m)expoS,(-imO I (/^"(O)'^' dl^fe) ■ ■ ■ dl?(e. 



Then 
with 



' ( — 77, 7r] 



W^ll^ = , (49) 

f^'hO = (expoS,(ifc7.0) ° (/^nO)'^' X l^^^^iO , 

where ^ = (^i, and /^^O = /(^i) ■ ■ ■ /(e,)- 

The two following results provide the asymptotic behavior of each term of the sum 
in (H7I) and of Wj^if\ respectively. They are proved in Sections l6.ll and W2\ respectively. 
The first result concerns the terms with long memory, that is, with q < Qc- The second 
result concerns the terms with short memory for which q> q^. 

Proposition 6.1. Suppose that g G {1, . . . , qc\- Then, as j — ?■ oo, 

(77(^(«)+^)w^t,., m,kez)^ {{nO)y^' Yif\ m^kez) , (50) 
where Y^'j^^ is given by [3l\ ). 
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Proposition 6.2. We have, for any A; G Z, as j oo, 

= Op{lf) ■ (51) 

It follows from Proposition 16.11 that the dominating term in ( l47l) is given by the chaos 
of order q = go- Now, since d{qo) > by ( l32i) . we get from Proposition 16.21 that, for all 
{k, m), as j — )■ oo, 

This concludes the proof of Theorem 15.11 

6. 1 . Proof of Proposition \6.1\ 

We first express the distribution of {w'^''^^^^, m, A; G Z} as a finite sum of stochastic 
integrals and then show that each integral converges in L^(fi). 

Lemma 6.1. Let g G N*. For any j 

[g/2] 

E (52) 

s=-[q/2] 

where [a] denotes the integer part of a, and for any g G N*, s G Z, 



w^Sr^= / W(77'C)/m,fc(C;j,g)dw^(Ci)---di?(g, (53) 

where fm,k{C'J,Q) is defined by (setting ^ = --fJ^C) 

Mfcl7j^,J,gj-7,- {l-expoS,(-iO}^ ^^'^ ■ ^ ' 

and where 

r^''-") = 1^ ^ (-^' -TT + 2s7r < ^ < TT + 2s7r| . (55) 
Proof. Using ( l49l) . with j replaced by j + m, and ([9]), we get 

By fl5^ . we thus get 

W^im,k = [ iff^Al^i: J, q) dW?(ei) ■ ■ ■ ^W{Q (56) 

>''£6f-7r,7rl9 



(fldi) 



£6(-7r,7r]' 

/„,fc(C;j, g) dW?(Ci)---diy(g, 

Ce(-73'f>73'r]'' 
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where we set C = 'jj^ (see Theorem 4.4 in [l7|). Observe that for all ( E (—7^71,7^7?]'', 



-T^lj - 2[g/2]7r7j- < -q-fj7i < J^O < Qlj-^ < -^Ij + '^[q/'A'^lj ■ 



i=l 



The result follows by using that for any C, G (— 7j7r, 7j7r]'', there is a unique s = — [g/2], . . . , [q/2] 
such that C/7i e T^^'"). □ 

Proof of Proposition In view of Lemma we shall look at the L^(fi) convergence 
of the normalized W^^'^'*^ at each value of s. Proposition 16.11 will follow from the following 
convergence results, valid for all fixed m, G Z as j — t- 00. For s = 0, 



1: 



(57) 



whereas for other values of s, namely for all s E { — [q/2], . . . , —1, 1, . . . , [fj/S]}, 



where d{q) is defined in (I3T1) . 

We now prove these convergence using the representation fl53|l . By ([1]) and |1 — e^'*'| > 
2|A|/7r on A G (— 7r,7r), we have that 



/(A)<(-) ||rlloo|Ar^^ AG[-7r,7r] 



(59) 



By definition of T^'?''') in ([55]), we have, for all C G 7jr(«'^), 7^"^ Xli -27rs G (-vr, tt]. Hence 
using the (27r)-periodicity of /ij+m, we can use @ for bounding /ij+m(7j~^ X]j Ci)- With 
the change of variables ( = 'jj^ and (|59|1 . for all C G 7jT('^''*) and j large enough so that 



(60) 



where Cq is a positive constant and 



9{C,t)= 1 + 



The squared L^-norm of g{-;t) reads 



EC- 



i=l 



-a-K 



nic. 

j=i 



EC. 

i=l 



-2a-2K 



4=1 



j = l 
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We now show that Lemma 18.41 apphes with Mi = 2a + 2K, M2 = and (3i = 2d for 
i = 1, . . . ,g. Indeed, we have M2 — Mi = —2a — 2K < —2a < —1. Further, for all 
£ = l,...,g — 1, we have, by the assumption on d, 

J]A = 2rf(l + g-£) > {1 + q - e){l - 1/q) = q - e + {£ - l)/q > q - e . 

Finally, since a > 1/2, one has M2 - Mi + q = -2a - 2K + q < q - 1 < f3i. 

Applying Lemma we get J(t) — )■ as |t| — )■ 00 and J(0) < 00. Thus, if s 7^ 0, one 
has t = 2n'jjS — )■ 00 as j — )■ 00 and hence we obtain (!58|l . If s = 0, then t = 2'K'^jS = 
and using the bound (!60l) . J(0) < 00, and the dominated convergence theorem, we have 
that the convergence (|571) follows from the convergence at a.e. ^ G M'^ of the left hand side 
of (|60l) . which we now establish. Recall that fm,k is defined in (J54l) . By ([6]), (II]) and the 
continuity of /* at the origin, we have, as j — )■ 00, 



and for every £ = !,■■■ ,q 



^ lU'h^{-f^{Ci + ■■■ + (,)) , 



if'fiCih,) = 77'1i - e-'^'/^^r''/*(0/7,) ^ /*(0)|0l 

Hence 77^'^^^^^^Vm,fc(C; J, 9, 0)lr(.,s)(77^C) converges to 

^ (i(Ci + --- + C.))^ICi|^---|C.I'^ ■ 

This concludes the proof. □ 

6. 2. Proof of Proposition \6.2\ 

We now consider the short-range dependence part of the wavelet coefficients {Wj^k) 
defined by (l29ll and (HHll . These wavelet coefficients can be equivalently defined as 



-2d 



wjP^Y.lf'i-^k-DA-'Yr, (61) 

where we have set 



g>gc+i ^ 

Using Lemma [4.11 since fl26|) holds and {Hq{Xe)}e^z are uncorrelated weakly stationary 
processes, the process {A^Y^^^''}g(zz is weakly stationary with spectral density 

/"^HA)= E %(-HA), AG(-vr,vr). 
<?>gc+i ^' 
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By Lemma Sm^ii), we have that ||/^*'f''"+-^^^||oo < oo. Using that \\gi-kg2\\oo < ||5'i||oo||5'2||i 
and II /111 = 1 by assumption, an induction yields 



SUp||/(*^)|U<||/(*^^^+'»||oo 

g>qc 



Hence, by ( 126|) . we get ||/*-"^^^||oo < oo. It follows that, for Wj^^^ defined in (|6T]) . there is 
a positive constant C such that. 



m^P"] < ll/^'^^lloo / \hf\X)\'dX < C / lAI-^^ \%{X)\\\X = 0(7f ) 



where we used dl]) with M > K and a > 1/2. This last relation implies f lSTj) and concludes 
the proof of Proposition 16.21 □ 



7. Proof of Lemmas 15.11 and 15.21 

7. 1 . Proof of Lemma I5.il 

Let us first prove that if ^ e then Z^JV) exists. Indeed, by Definition Z^^f{e) 
exists if ^ 

- — — -^j^ -^dui ■ ■ ■ dwq < oo . (62) 

1^1 + ■ ■ ■ + Mqp-" |Mi ■ ■ ■ Mgl 

Use now Lemma \K3\ with /3i = ■ ■ ■ = /9g = —2d and f{x) = \6{x)\'^ /\x\'^^ and deduce that 
Condition f l62l) is equivalent to 



where 



r / |^(s)ns|''-i*2^ds<oo. 



(63) 



Note that the conditions < d < 1/2 and < q < 1/(1 — 2d) ensure that F is finite. 

g,d ■ 

^ 



Further, Relation ( 1^ implies 6 G S^^\ 

We now prove that for any m, /c, hoo,m,k £ <Sjf2 when K G {0, ...,M}. By Defini- 
tion (HOD of /loo,m,fc 

Hence 

/ \hoo,rnA^)\'\sr'~'''-'''ds = J^ [ \h^{-f^s)\'\sr'-"^'-'''ds 

Jw Jr 
Set V = and deduce that hoo,m,k £ '5^'^'* is equivalent to 

JR 
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Assumption (fT2!) implies that 



\2M 



il + \v\) 



2M+2a I 



\q-l-2qd-2K 



Since M > K and q{l-2d) e (0, 1) then 2M+q-l-2qd-2K = {2M-2K)+q{l-2d)-l > 
-1. Further a > 1/2 and g(l -2rf) G (0, 1) imply that 2M-2M-2a + (g-l-2grf-2ir) = 
-2a -2K + q{l - 2d) - 1 < -1. Then 



\hr 



|2| ,,-l-2,d-2if^^ < OO . 



holds and hoo,m,k e Sj^^J . 



7. 2. Proof of Lemma \5.2\ 

Let at{ui, ■ ■ ■ ,Uq) denote the kernel of the integral in fj42|) defining Z^^^^ and suppose 
we can exchange the order of integration and write 



Then condition (H3l) gives 



M.1 



at{ui, ■ ■ ■ ,Ug)d{t)dt 



dW{ui)---dW{ug) . (64) 



K-l 



(it(Mi + ■ ■ ■ + Ug)Y 



e{t) dt = / e"("^+-+"^) 9{t) dt = eo Eg{u) 



showing that ([MD equals zfJiO) defined in dnHD- It remains to justify the change of 



order of integration in (1641) by using a stochastic Fubini theorem, (see for instance 20 
Theorem 2.1]). A sufficient condition is 



, ■ ■ ■ , Ug)dui ■ ■ ■ duq) ^^"^ dt < OO 



This condition is satisfied, because setting v = tu, we have 

2 



£=0 



|i(Ml + --- + M,)|-2^|Mi---Mg|-^^ d% 



-2d 



<|^|2X+2d-, f (l+|^^ + ... + ^J)-2X|^^...^^|-2d^,^_ 
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8. Auxiliary lemmas 

The following lemma provides a bound for the convolution of two functions exploding 
at the origin and decaying polynomially at infinity. 

Lemma 8.1. Let a > 1 and /3i,/32 € [0, 1) such that /3i + /32 < 1? (^nd set 

g,it) = \t\~^^{l + \t\f^-" . 

Then 

sup((1 + |m|)° gi{u-t)g2{t)dt] < oo . (65) 

Proof. We first show that 

J{u)= [ gi{u-t)g2{t)&t= I \u-t\-^^{l + \u-t\f^-''\t\-^^{l + \t\Y^~''At 
Jr Jr 

is uniformly bounded on R. Using the assumptions on /3i,/32, there exist p > 1 such that 
/3i<l/j9<l — /32- Let q be such that l/p + l/q = 1. The Holder inequality implies that 

j(m)p«< [ |t|-p^i(i + |t|)p^i-p"dt X [ |t|-''^^(l + itD^^^-^^dt . 

Jr Jr 

The condition on a, f3i, I32,p and the definition of q imply that these two integrals are finite. 
Hence sup„ J{u) < oo. 

We now determine how fast J{u) tends to as m — > oo. Observe that, if |t — < \u\/2, 
then \t\ > \u\/2. By splitting the integral in two integrals on the domains |t — < \u\/2 
and |t — m| > \u\/2, we get J{u) < Ji{u) + J2(m) with 

Ji(n) < (|n|/2)-^^(l + |m|/2)^2-" f \u- t\-^'{l + \t- u\f'~''dt , 

Jr 

and 

J2{u) < {\u\l2)-^^{l + \u\l2f^-^ I \t\-^\l + \t\Y^-^dt . 

Jr 

Now, as \u\ — 7- oo, we have Ji{u) = 0(|m|~") for i = 1,2, which achieves the proof. □ 

The next lemma describes the convolutions of two periodic functions that explode at 
the origin as a power. A different definition of convolution is involved here (see ( l33l) ). 

Lemma 8.2. Let (/3i,/32) £ (0,1)^- Let gi, g2 be (27c) -periodic functions such that gi{X) = 
|A|^^* (7*(A), i = 1,2. Each g*{X) is a (27^) -periodic non-negative function, bounded on 
(— 7r,7r) and positive at the origin, where it is also continuous. Let g = gi g2 o-s defined 
m Then, 

• If f3i -\- 132 < 1 , g is bounded and continuous on {—11,71), and satisfies g{0) > 0. 
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. ///3i+/32>l, 

^(A) = \\\-^^^+^--'^g*iX) , 

where g*{X) is bounded on (— 7r,7r) and converges to a positive constant as A — ?■ 0. // 
moreover for some /3 G (0, 2] such that /3 < /3i + — 1 and some L > 0, one has for 
any i G {1, 2} 

k*(A)-(7*(0)|<L|A|^VAG(-7r,7r), (66) 
then there exists some L' > depending only on L, Pi, such that 

|^7*(A)-(7*(0)| <L'|A|^ VAG(-7r,7r). 



Proof. By (1331) and (27r)-periodicity, we may write 

g{\) = r gi{u)g2{\ -u)du= T |{A - u}\-^'gl{\ - u) \u\-^'gl{u) . (67) 

J —n J —IT 

Let us first consider the case Pi + (32 < 1. We clearly have (^(0) > 0. To prove that 
g is bounded, we proceed as in the case of convolutions of non-periodic functions (see the 
proof of Lemma [HTT]) . namely, for p, q such that /3i < 1/p < 1 — and 1/p + 1/q = 1, the 
Holder inequahty gives that 

Ikfo? < 11^211^ < \\9l\L \\9;\\lo r Itr^'dt X r Itr^Hlt < oo . (68) 

J —TT J —n 

For any e > and i = 1,2, let g^^i be the (27r)-periodic function such that for all A G (— vr, it), 
ge,i{^) = l(-,,c)(A) gi{\) and let g,^i = g-i-g^^i. Then g = ge,i-^ge,2 + 9e,i'^9e,2 + 9e,i-^ 9e,2 + 
9t,i 9e,2- Since g^^i is bounded for i = 1,2, we have that (7^,1 * ge,2 is continuous. On the 
other hand, using the Holder inequality as in f l68|) . we get that H^/e,! * (7e,2||oo5 ||^e,i*5'e,2||oo) 
*^7e,2||oo tcud to zcro as e — 7- 0. Hence g is continuous as well. 
We now consider the case (3i + (32 > 1. Setting v = u/X in ([HZD, we get, for any 
AG [-7r,7r]\{0}, 

^*(A) = \Xf^+f^--'g{X) = [ l(_./|,|,./|,|)(t;) \{{l-v)}^\-^^\v\-^^g*i{X{l-vM{Xv) dv , 

where for any real number x and A 7^ 0, {x}x denotes the unique element of [— 7r/|A|, 7r/|A|] 
such that X — {x}x G Z. Take now |A| small enough so that 7r/|A| > 2. Then, for 
any v G (— 7r/|A| + l,7r/|A|], we have |{(1 — v)}x\ = |1 — f | > |1 — |f|| and, for any 
V G (— 7r/|A|, — 7r/|A| + 1], we have 

\{(l-v)}x\ = |l-t;-27r/|A|| = 27r/|A| +t; - 1 > -t' - 1 = |1 - . (69) 

Thus we have l(_jr/|A|,7r/|A|)(^^) |{(1 ~'y)}A|~^^| < |1 — l"^!!"^^ ^oi all t> G M. We conclude 
that for |A| small enough, the integrand in the last display is bounded from above by 
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|1 — It'll '^^If I ^'^ Iki ||oo||fi'2lloo, which is integrable on t> G M. Hence g* is bounded, and by 
dominated convergence, as A — )■ 0, 



g-iX) ^ gmgm / |1 - > 



(70) 



We set g*{0) equal to this hmit. 

Suppose moreover that gl,g2 satisfy ( l66l) . We take gKO) = g2{0) = 1 without loss of 
generality and denote rj(A) = |5'*(A) — 1| for i = 1,2. Then r(A) = |5'*(A) — 5'*(0)|, where 
g*{0) is defined as the limit in ( !70|) , is at most 



|1W|AU/|A|)(^) m-v)}xr^^\v\-^-gliX{l-vM{Xv)-\l-v\-^^ \vr^-\ dv. 
Setting ^*(A) = (^*(A) - 1) + 1, we have r<A + Bi + B2 + C with 

A{X)= [ |{(l-,;)},|-/^>|-/^2-|l-^;|-/5l \v\-P^\ dv , 

B,{X)= [ l^_^/ixi,./\xi){v)\{il-v)}x\-^^v\'^'n{Xv)dv , 
where is (1,2) or (2,1), and 

C(A)= / l(-./|A|,./|A|)(^) \{il-v)}x\-^^\v\-^'n{Xil~v))r2{Xv)dv . 

JR 

Since {(1 — f )}a = 1 — for v G [— 7r/|A| + 1, 7r/|A|) and A large enough, we have 

A{X) = I |l-w|-^^ \v\-^^ dv 

J{-n/\\\,^/\\\Y 

I'-tt/\\\ + 1 

+ / \\{{l-v)}xl^'\v\-^' -\l-v\-^'\v\-^''\ dv . 

J-7t/\\\ 

The first integral is 0(|A|'^^"'''^^~^). Using fl69l) . the second line of the last display is less 
than 

/ [\l-v\-^'\v\-'^'' + \l + v\-'^' v-f^'] di; = 0(|A|^^+^^) . 

Jtt/\\\-1 

We conclude that as A ^ 0, A{X) = 0(|A|^i+^2-i)_ Moreover using that ri(A) < LlAI'' and 
/3i + /32 — /3 > 1, we have -Bj(A) = 0(|A|^) for i = 1,2. The same is true for C since ri and 
r2 are also bounded on M. This achieves the proof. □ 

Lemma 8.3. Let p be a positive integer and / : M — )■ IR+. Then, for any P G M^, 

/ fiyi + --- + y,)f\h\^^dy,---dy, = Tx [fis)\sr'+^'+-~^^''ds, (71) 
Jri -^j^ Jr 
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where, for all i & {1, ■ ■ ■ , q}, Bi = f3i + ■ ■ ■ + f3q and 

g 



i=2 



\t 



q-i+Bi 



1 1 . 



(We note that V may he infinite in which case (fTiy holds with the convention oo x = 0^ 



Proof. Relation (I7T1) is obtained by using the following two successive change of variables 
followed by an application of the Fubini Theorem. Setting, for all i = 1, • • • , g, Ui = 
Yl'j=iyjy get that Hi = Ui — Wj+i for i < g and Uq = Uq. Then the integral in the 
left-hand side of (1711) reads 





q-1 




i=l 



dui ■ ■ ■ duq . (72) 
The second change of variables consists in setting, for alH = 1, ■ ■ ■ , g, = 11}=! ^i- Then 



dui---duq = I n*r' 1 dtf-dl 
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q-l 

n 



q-l 



^q-l 



Ui - Ui 



i=l i=l \i=l 

and l^gl = riLi \^i\^''i so ^^^^ (1721) becomes 

/ fiti) n nil - • • ■ dt 

i=l i=2 

which by Fubini Theorem yields the required result. 



-t. 



vi=2 



□ 



Lemma 8.4. Let a G M and q he a positive integer. Let /3 = ■ ■■ G (— oo, 1)'^, 
Ml > and Ma > -1 s^ic/i i/iat Ma - Mi < -1. ^sswme t/iai g + Ms - Mi < YlUi A, 



and i/iai /or any £ G {1, ■ ■ ■ , g — 1}, f3i > q — i. Set for any a G M, 



Jqia; Ml, M2; 13) 



|S,(C)-a 



Ah 



;i + is,(c)-ai)^^^nic. 

i=l 



Then one has 



sup(l + |a|)^-''+^-i^'Jg(a;Mi,M2;/3) < 00 

aeK 



(73) 
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In particular, 

J,{0; M,, M2; < 00, 

and 

Jg{a; Ml, M2; (3) = 0(|a|-(i-^+S-i'^')) as a -> 00 . 

Proof. Since Jq{a; Mi, M2; ■ ■ ■ , (3g) = Jq{—a; Mi, M2; (3), we may suppose a > 0. By 
Lemma 18.31 



J,(a; Ml, M2; /3i, ■■ ■,/?,) = r / ^ — ds 

Jr U + p — c^U 



where 

i=2 



The conditions on /3j's, Mi and M2 imply Jg(a; Mi, M2; /3i, ■ ■ ■ , < 00 for all a. To 
obtain the sup on a > 0, we set v = s/a. Then, denoting S = XlLi Z^*' 

J,(a;Mi,M2;/3) =Ca«+^'2-^ / |t; - + a|t; - iD-^^j-^+^^-^^dt; , (74) 

Jr 

where C is a positive constant. We separate the integration domain in two. Suppose first 
that |f — 1| < a~^. Then in this case we have (1 + a\v — < 1. Since \v\ is bounded 

on the interval |f — 1| < a^^ for a large then as a — )• 00, 

f |t;-l|*^2(i+^|^_i|)-Mi|^|-5+(,-i)^^ = o( [ \v- ll'^'dv) = 0(a-^-^2) , 

J|i)-l|<a-i \J\v-l\<a-^ J 

Now suppose that \v — 1 > a~^\. Then (1 + a\v — < {a\v — and 



V 



(ii>2 1^ - ir^-*^>r^+(-^)dt; + /i/2<M<2,i.-ii>.-. 1^ - ir-''^\v\-'^^^-'^dv 

The first integral concentrates around v = 00, the second around f = 1 and the third 
around f = 0. The first integral is bounded, the second is 

0{j \v- l\''^-^'^dv) = 0{a^'-^'''-'), as a ^ 00, 

and the third is bounded. Therefore we get 
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since M2 - Mi < -1. Thus gives 

Jg(a;Mi,M2;/3) = 0(a~^+«"^) as a ^ 00, 

yielding the bound ( j73l) . □ 
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